DILUTION OF COSMOLOGICAL DENSITIES 
BY SAXINO DECAY 

David H. Lyth* 

School of Physics and Materials, 
Lancaster University, 
Lancaster LAI 4YB, 
U.K. 

May 1993 



Abstract 

Saxino decay can generate significant cosmological entropy, and hence dilute theo- 
retical estimates of the present mass density of a given particle species. The dilution 
factor depends on the saxino and axion masses, and is constrained by the requirement 
that saxino decay should not affect nucleosynthesis, as well as by the usual requirement 
that the axion density be less than the critical density. The latter constraint is eval- 
uated carefully, under both the Harari-Sikivie and Davis proposals about the emission 
spectrum from axionic strings. Uncertainties are carefully evaluated, points of principle 
are addressed, and with an eye to future numerical simulation the spacing and typical 
oscillation wavelength of the strings are represented by parameters varying in the range 
1 to 3. Within the constraints, the entropy dilution varies from 1 to 10^"'. Only saxinos 
originating from thermal equilibrium are considered, so that more dilution might arise 
from non-thermal saxinos. 
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1 Introduction 



According to current ideas, the present mass density of the universe is practicahy equal 
to 1 (in units of the critical density), baryons contributing ~ .04 to .07, and one or 
more species of dark matter particle making up the rest. Given a suitable model of the 
fundamental interactions one can calculate the present density Ox of a given stable particle 
species. The calculation usually proceeds in two stages fll. First one calculates the particle 
number density at some initial epoch, defined by its temperature Tj, after which the particle 
number is conserved. Then the number density is evolved forward to the present epoch when 
the (photon) temperature is Tq = 2.74 K. 

As long as thermal equilibrium is maintained, the forward evolution can be performed 
using entropy conservation. Entropy can however be generated by non-equilibrium pro- 
cesses, and this lowers the prediction for Qx- A particular example is the decay of a 
particle species which dominates the energy density of the universe during some era prior 
to its decay. The Standard Model does not contain such a particle species, but its as yet 
unknown extension might. 

Any acceptable extension must include a natural explanation of CP conservation by the 
strong interaction. Although it has recently been seen to be rather vulnerable to Planck 
scale corrections ^, the most promising explanation is still a spontaneously broken global 
symmetry of the type proposed by Peccei and Quinn [Q, ^] . Along with this symmetry comes 
a Goldstone boson, the axion |5|, [l|, The axion is practically stable and has extremely 
weak interactions, constituting cold dark matter. From accelerator physics and astrophysics 
its mass has an upper bound ma < 10^^ eV, and from the cosmological requirement < 1 
its mass has a lower bound ma > m^i^. The lower bound mmm is difficult to calculate, but 
with no entropy production it is probably not far below 10~^ eV. In that case there is only 
a narrow allowed window for ma, and the axion must make up a significant fraction of the 
dark matter. 

Despite the fact that no supersymmetric particle has yet been observed, it is fair to 
say that an increasing number of particle theorists feel that an acceptable extension of 
the Standard Model should also respect (low energy) supersymmetry Q. One reason for 
this feeling is the successful prediction of sin^ 9w by supersymmetric grand unified models, 
and another is the fact that the rival technicolour models are coming under pressure from 
ever more accurate measurements, notably at CERN. Supersymmetry requires g that the 
known particle species have as yet undiscovered superpartners, with masses of order .1 to 
1 TeV. It also requires that the axion be accompanied by a spin 1 /2 partner called the 
axino, and a spin partner called the saxino [10, 11, [l^. The axino might have the typical 
mass mentioned above, in which case it would decay (perhaps with significant entropy 
production), or it might be much lighter in which case it could constitute 'warm' dark 
matter. (More complicated possibilities arise if there is a very light gravitino [13|, which 
are discounted in the present paper.) The saxino, on the other hand, is definitely expected 
to have a mass in the usual range, 100 GeV < msax ^ 1 TeV. 

Kim 1 12] has pointed out that saxino decay could generate significant entropy, and in the 
present paper the amount of entropy production is carefully calculated as a function of the 
axion mass ma, the saxino mass msax and the maximum temperature after inflation ('reheat' 
temperature) Tj-eh- Then, the cosmologically allowed region of this three parameter space is 
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delineated, by the requirement that saxino decay must not interfere with nucleosynthesis, 
as well as by the requirement that the axion density satisfies ila < 1- As described in the 
conclusion, the results have a number of cosmological implications. 

Throughout the units are h = c = 1, mpi = 1.2 x 10^^ GeV is the Planck mass, a{t) 
is the scale factor of the universe at time t, H = a/a is the Hubble parameter, and Hq = 
100 km s^^ Mpc^^/i. is its present value. 



2 Entropy and particle densities 

For reference we recall some basic facts about radiation (relativistic particles) in the early 
universe ||^] . According to the Standard Model the radiation is in thermal equilibrium when 
the temperature T exceeds 1 MeV, with practically zero chemical potential. Accordingly its 
density is 

where is the effective number of relativistic particle species Its entropy density s is 

Apr 2^2 

With three massless neutrino species, 5* = 10.75 for IMeV <T< 100 MeV, g^ = 17.25 for 
100 MeV < T < Aqcd and 5* = 61.75 for Aqcd < T < 2GeV where Aqcd ^ 200 MeV. 
At higher temperatures the Standard Model predicts a leveling out to 51* ~ 100 above T ~ 
100 GeV. However, extensions of the Standard Model have additional particles, leading to 
a bigger value of 5*, and in the minimal supersymmetric extension 5* levels out at 5* = 229. 
Extending the Standard Model does not usually have any other effect on the formulas that 
we are reviewing in this section, except for the possibility of entropy production that is our 
central concern. 

In thermal equilibrium the entropy S oc a^s in a comoving volume is constant. This 
implies that 

a^g^{T)T'^ = constant (3) 

After the neutrinos decouple at T ~ 1 MeV they are not in thermal equilibrium but their 
entropy continues to be conserved, and the same is true of the photons when they too 
decouple at a much later epoch. Thus, according to the Standard Model the total entropy 
in a comoving volume is conserved, and a standard calculation |Q| gives the present entropy 
density sq = 2938 cm~'^. This allows one to calculate the present mass density of any 
particle species X, whose number density is conserved after some epoch Tj > 1 MeV, given 
the number density at that epoch. Indeed, in units of the present critical density SHqItipi/Stt 
it is obviously given by 

_ Sir msQ rij 

= 2.Uh-^^^ (5) 
10 eV Si ^ ' 

This relation has been used to calculate the baryon density VLb and the density J7x of various 
dark matter candidates, in terms of theoretically calculated initial number densities. 
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The above analysis may fail if there is a non-relativistic particle species with a lifetime 
in the range tj < r < to- However, the required modification is very simple provided that 
the decay products thermalise quickly on the Hubble timescale. Then the radiation density 
and entropy density are still given by Eqs. (|T|) and (^), and the only effect of the decay is 
to reduce Ox by a factor A = So/Si, the increase in the entropy of a comoving volume 

Neither the Standard Model or its minimal supersymmetric extension contains such a par- 
ticle species, but for extensions which include Peccei-Quinn symmetry the saxino can be 
such a species, as we now discuss. 



3 Entropy production from saxino decay 

The supersymmetric generalisation of the usual axion-gluon interaction is ll 



C = -^^W,W' (7) 

where the superfields are <I> corresponding to the axion and saxino (spin 0) and the axino 
(spin 1/2), and Wi corresponding to the gluon (spin 1) and the gluino (spin 1/2). This 
reproduces the usual axion-gluon interaction, and also gives a saxino-gluon interaction and 
other interactions, 

^ = ^ (V'aG^.G'^" + ^saxG^.G^'^ + . . .) (8) 

The saxino is kept in thermal equilibrium by reactions like qq sg and gg ^ sg above a 
temperature |12| 

r.ec„„p = 10"Oev(^^)'(^)"' 0) 

At an energy scale of order 10^^ GeV, ac — 1/20 in the minimalQ supersymmetric standard 
model 1 14], so that 



f: 
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Tdecoup =.8x10^^ GeV [j^^^ ) (10) 

After decoupling the saxino number density rigax is given as a fraction of the entropy density 
by 

r = ^ = = 1.2 X 10"^ (11) 

s 3.65* 

(setting g^, = 229 corresponding to the minimal supersymmetric standard model). 

In order to achieve thermal equilibrium the reheat temperature must satisfy T^eh > 
Tdecoup- Unfortunately, even the order of magnitude of Treh is unknown. Bounds on the 
cmb anisotropy, as well as the COBE detection, imply that it is less than 10^^ GeV [15, 16|. 
If there is a gravitino with mass rUg < 1 TeV nucleosynthesis implies that it is less than 

^ Since Peccei-Quinn symmetry is assumed the extension cannot in fact be minimal, but this estimate of 
Qc and the estimates of g* used later will hopefully still be adequate. 
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10^^ GeV, but that bound rapidly goes away as nig is increased and provides no additional 
constraint if rUg ~ lOTeV [0]. 

In what follows we leave free parameter, and assume that if Tj-ch < Tdccoup 

the saxino density is negligible, considering only saxinos of thermal origin. A significant 
non-thermal saxino density might originate as an inflationary fluctuation or be radiated 
from axionic strings, possibilities which should certainly be investigated. 

Continuing the story of the thermal saxinos, they become non-relativistic at a tempera- 
ture T ~ msax) and if they live long enough they dominate the energy density of the universe 
below the temperature 

Tsaxeq = ^rmsax = 1-6 GeV(^^^) (12) 

(The subscript 'saxeq' indicates that the saxino density is equal to that of the radiation at 
this temperature.) The dominant decay mode of the saxino is into two gluons, with lifetime 
r-i given by |l| 

9 

-p ^c"^sax (13) 



128^3/^2 

and saxino domination will occur for a significant period if tsaxeq ^ r~^. 

With this assumption, the evolution of the matter and radiation densities can easily 
be worked out. Approximate analytic results are available on the assumption that 17* is 



T— — Tjgi^g^^Tgaxeq (15) 



constant [18, |l^, ^ |2l|, y, which are now briefly recalled. The temperature corresponding 

to tdecay is 

Tdeeay = .55g7^^'^ V^rup (14) 

and the saxino dominates the energy density until this epoch. The density of the 'new' 
radiation from saxino decay becomes equal to that of the 'old' radiation at a temperature 
r=, given by 

rr5 _ r 

decay -'saxeq 

Before this epoch there is no significant entropy generation, but after it the entropy in a 
comoving volume is proportional to . This continues until the epoch Tdecay; after which 
the saxino density becomes negligible and entropy production stops. Thus, if Tj < T= the 
entropy generation factor in Eq. (^ is 

A = [t^^ m < r=) (16) 

\ ^ decay / 

It depends strongly on Tj because only part of the entropy is generated after Tj. On the 
other hand, if Tj > T=, all of the entropy is generated after Tj, and A is given by 

Atotai = ' = ^ (T, > r=) (17) 

y ^ decay j J- decay 

The last equality follows from Eq. (^). It implies incidentally that the new radiation has 
caused the ratio Psax/Pr to fall to a value ~ 1, just before it falls off exponentially at the 

epoch Tdecay 
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These results are valid if g^, is constant during the interval Tjecay <T < Tgaxeq- Setting 
5* = 10.75, appropriate to the range IMeV <T < 100 MeV, gives 

r.„., = 52Mev(^)"(^)"^ (18) 

The three temperatures Tdecay ^saxcq srs plotted in Figure 1. If the saxino 

is heavy, Tjecay can exceed the temperature T = Aqcd at which g^: rises sharply to 62, 
but using this value would not change the results much. We have not investigated the 
intermediate case T^ecay < ^qcd < T'saxeqj but assume that the results above are still valid 
to sufficient accuracy. 

Substituting Eq. ( |l8|) into Eq. ( [l7|) gives the maximum entropy generation factor which 
can be supplied by saxino decay, 

A«>. = 3l(^)(^)"'' (19) 

The saxino must decay before it has a significant effect on nucleosynthesis. According 



to Scherrer and Turner |^ a safe bound is t < 10 s, which is equivalent to Tdocay ^ .3 MeV, 
or 

The region forbidden by this constraint is indicated in Figure 1. The corresponding con- 
straint on the entropy generation factor is 

AtotaK 5.2x10=^1^^^) (21) 

A further constraint can arise from the requirement < 1; which is our main concern for 
the rest of the paper, as is now investigated. 

4 The axion density 

The axion density has been the subject of numerous investigations. The situation as it was 



understood in 1989 is described in [PJ, Ml) and further work has been done since then [22, 



23, U, 25, 27, 1^. There are essentially two possibilities regarding axion cosmology. 
One is that the Peccei-Quinn field has been in its vacuum (ie., that PQ symmetry has been 
spontaneously broken) ever since the observable universe left the horizon during inflation. 
In that case axion cosmology is quite complicated as will be briefly recalled in Section 4.3. 

More likely is the opposite case, where Peccei-Quinn symmetry is restored in the early 
universe, through either finite temperature effects (after inflation) or the quantum fluc- 
tuation (during inflation). In that case strings form at the epoch when the symmetry is 
spontaneously broken, and the strings oscillate until domain walls form when the axion 
mass switches on at T ~ 1 GeV. If the string network 'scales' as one expects, can in 
principle be calculated uniquely in terms of /a (up to the entropy generation factor A^^), at 
least if the strings radiate nothing but axions. To do a proper calculation requires numerical 
simulations which are perhaps tractable while the strings are oscillating but which would be 
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very difficult after domain walls form. In the absence of such simulations, ila has been esti- 
mated in the literature on the basis of simplifying assumptions. The first of these, used by 
essentially all investigators, is that axion number is conserved after the domain walls form, 
so that Eq. gives f^a in terms of the axion number density just before wall formation. 
The second assumption is a statement about the shape of the spectrum of axions emitted 
by the strings. Two alternative proposals exist, one by Davis [29| and one by Harari and 
Sikivie [30|. In what follows these two proposals are critically assessed, and in each case J^a 
is carefully estimated. The treatment represents an advance on previous ones (except for 
|26| which is discussed below) in two respects. First, the spectrum of all axions present at 
a given epoch is calculated (not just that of those being emitted), which enables important 
points of principle to be addressed. Second, the result is given as a function of the string 
spacing, as well as (for Davis' proposal) the string oscillation frequency. Correspondingly 
there appear in the result two parameters 7 and /?, which are expected to lie in the range 1 
to a few, and which should be calculable in the forseeable future from numerical simulations. 
Pending such a calculation, we estimate Oa by allowing them to vary between 1 and 3. 



4.1 Axion cosmology 

First we need some basic properties of the axion It is a practically stable particle, 

which is massless for T ^ Aqcd — 200 MeV. At lower temperatures its interaction with 
the gluon leads to an effective potential 

V = (79 MeV)^(l - cos 6) (22) 

where 6 = V'a/Za and -^a is the canonically normalised axion field. (We ignore a possible 
generalisation of the above equations, 6 — > N9 and /a fa./^ with N an integer, because it 
probably leads to cosmologically unacceptable domain walls.) Small oscillations around the 
minimum of this potential correspond to free axions with mass ma given by /a?Tta = 79 MeV, 
or 

6.2(^-4-)" (23) 



10-6 eV ■ VlQi^GeV, 

Accelerator physics and astrophysics provide a lower limit rria < 10^^ eV, corresponding to 
/a>10io GeV m. 



The case of a homogeneous axion field 

As already mentioned, axion cosmology depends crucially on whether or not there are 
axionic strings. Before turning to the former case which is our main concern, we focus on 
the simplest possible example of the latter case. Namely, we assume that the axion field is 
spatial homogeneous. In this case an accurate estimate of Oa is provided by the calculation 
of Turner [ |32| , which we briefiy recall. It will form the basis of our later estimates for the 
more complicated string case. 

The axion mass switches on gradually as the epoch T = Aqcdj becoming significant at 
the epoch i defined by 

ma(t) = (24) 
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Before this epoch, the homogeneous axion field is time independent, and afterwards it 
oscillates around = 0. If |0| <C vr (a tilde on any quantity will always denote its value at the 
epoch t), the oscillations are harmonic with angular frequency rh and amplitude proportional 
to a~^/^ma corresponding to the presence of zero momentum non- interacting axions with 
number density 

n = ^mfje^ (25) 



(Note that n oc corresponding to conserved axion number ||3^.) The present axion 
density is therefore given by Eq. (^) once T is known. While the axion mass is switching 
on it is given by 

m^{T)/m^ = miiAQCB/Tf (26) 
Taking h = .5 ± .1 and Aqcd = (200 ± 50) MeV this leads to 



/a 
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and 

^a = 0.9xl0±-5(^^^^j A-'9' (28) 

The quoted uncertainty is essentially the one evaluate by Turner as arising from the uncer- 
tainties in (jr* and m^{T), slightly increased to allow for additional uncertainty arising from 
the values of h and Aqcd- (Bearing in mind the age of the universe we here discount the 
possibility /i ~ 1, which according to Eq. (^) would multiply by a factor .25.) To this 
level of accuracy, the expression is valid for \6\ < .Ovr, only the case |0| ~ vr requiring special 



treatment 27] 



The entropy production factor A is equal to 1 if Treh < Tdecoup or if Tgaxcq < Tdecay where 
the four temperatures are defined in the last section. Otherwise it is given by Eqs. ( p!6|) and 
([l7|), with Tj = T. From the left hand column of Figure 1, one sees that unless the saxino is 
rather heavy T is bigger than T=, so that the total entropy generation factor is experienced 
by the axions, given by Eq. (|T9|). 

Eqs. (H^) and (p8|) are derived under the standard assumption that T occurs during 
radiation domination. From Figure 1 one sees that if the saxino is rather heavy T can occur 
during (saxino) matter domination, in which case there are well defined correction factors 



|19|. In the present context these are not very significant, and will be ignored. 

The expression for given by Eqs. (|l6|), (p^) and (|2l) does not agree with that of 



Kim [12|. The reason is that he used formulae developed by Lazarides et al [19| in a 
different context, where the decaying object is not a saxino and where T falls in the interval 
Tdecay < T < T= throughout the relevant portion of parameter space. Figure 1 illustrates 
that this is far from true in the present case. 



4.2 The string scenario 

The above calculation provides a basis for calculating f^a on the assumption that there are 
axionic strings, if axion number is conserved after domain walls form (the validity of this 
assumption is briefly addressed at the end of the present section) . By looking at the typical 
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spatial and temporal variation of 9 it was argued in |^6[ that the walls form at about the 
temperature ~ T, and this estimate will be used in what follows. (More precisely it was 
argued that the formation temperature is lower by an insignificant factor 7"^/^ where 7 is 
the string spacing introduced below.) Substituting Eq. ( |25[ ) into Eq. ( |28| ) gives an expression 
for in terms of the number density n just before wall formation. 

1.18 



'a 



°-'"'°^ '(l0^4^) (2--a Va-^) (29) 



To estimate n, the simplest assumption is that 9 is as homogeneous as it can be, subject 
to the requirement that it changes by 27r as one goes around a string. This implies that the 
typical value of 9 is of order 1 radian, and to a rough approximation it also implies that 
the space and time dependence of 9 can be ignored in the equation of motion. After wall 
formation there is an a roughly homogeneous axion field oscillating with initial amplitude 
^ ~ 1, giving h ~ ^m-a/a ^'^d 

f \ 1-18 

Ja. \ A-l 



As we now discuss, this simplest estimate was queried by Davis [29|, who argued that axion 
radiation from strings plays a crucial role. 

Energy loss from the string network 

Assume that there is a network of axionic strings, with at least one piece of string passing 
through a typical Hubble volume. On scales bigger than the horizon the string network 
expands with the universe, and on smaller scales string annihilation occurs. We assume 
that the string network has the 'scaling' property, whereby it looks the same at each epoch 
when viewed on the Hubble scale. To be precise, we assume that the energy density of the 
string network can be written 

/Jst = ^Vt' (31) 

where /x is the string energy per unit length and 7 is time independent. Roughly speaking, 
7 is the string spacing in units of the Hubble distance H~^. 

The energy per unit length is dominated by the axion field around the string, rather 
than by the string core. The latter has thickness /a~^ and energy per unit length /a^, but a 
static axion field out to a distance R gives energy per unit length 27r/a^ ln(i?/a). As we shall 
discuss in a moment the axion field far away from the string is oscillating, but as a rough 
estimate we can set R equal to the string spacing {H'y)~^. Then n has only logarithmic 
time dependence, and at the relevant epoch T ~ 1 GeV it is given by ^ = 2'Kf^rj where 
T] = ln(/a/7i^). With /a ~ 10^^ GeV and 7 ~ 1 this gives the estimate r/ ~ 70, which we 
use in what follows. 



Prom Eq. (|3l| ) one can calculate the rate at which the string network loses energy |29|. As 
in the homogeneous case we assume radiation domination throughout the string oscillation 
era. One then has H = l/2t and Eq. (|3ll ) becomes 



Pst = '-I'k'vt--' (32) 
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The condition that comoving strings lose no energy is pst oc a^"^ . To maintain the scahng 
solution, energy conservation therefore requires that a comoving volume of the string 
network emits energy at a rate a^R where R is given by 

R=-a = -7 /ar/t (33) 

All of this energy is assumed to be gained by the axion field. In the absence of the string 
network this field would correspond to a collection of massless non-interacting axions with 
energy density pa oc a~^, so energy conservation requires 

a-^ = i^ = |7Va.r3 (34) 
Estimating the axion number density 

So far everything is rigorous, given the scaling assumption. To go further one has to 
introduce the concept of an 'axion' with well defined momentum and energy, as opposed 
to just an 'axion field'. Such an axion corresponds to field configuration which is a plane 
wave, and this clearly requires that the wavelength is much less than the string spacing, or 
equivalently that the wavenumber uj = k/a is much bigger than 

Wmin = = l-K-fH (35) 

a 

For a plane wave with amplitude 6, the energy density pa is related to the number 
density n by /Ja = i^?T-a where lo = k/a is the angular frequency. Let us define the spectrum 
T'a of the axion energy density as the contribution to it from unit interval of In/c, 

Pa= / mt,k)— (36) 

""=min ^ 

then 

f°° ^ , , N fa\ dk , 

"■mm ^ ' 

One can also define the spectrum of the axions which are being emitted during some 
small time interval dt, which we write as dtVemis- It is related by energy-momentum con- 
servation to the axion spectrum (cf. Eq. (p^), 

Vemisit, k) = a-^- (a^V^it, k)) (38) 



Integrating over k gives the energy conservation constraint Eq. (34), 

dk 



r VemUt,k)^ = ^l'U'vt-' (39) 

Jk ■ k 2 



Given a hypothesis about its shape, this equation determines Vemis, and integrating Eq. (^ 
then gives Va which gives the number density through Eq. (I 
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Two proposals have been made about the shape of the spectrum. According to Davis 
p9| , axion emission is caused by smooth oscillations of the string, with wavelength 

roughly f3~^ times the string spacing, where P is between 1 and a few. This corresponds to 

remis{t,k)=f{t)6{k-h{t)) (40) 

where 

iO* = — = PUJmin = —— (41) 

a t 



Inserting this expression into Eq. ( |38|) gives 

V^it,k) = 7r7]-f^UH-^e{k - k,{t)) (42) 

and then Eq. ( |37D gives 

duj 



n 



-Vl'fj-' r % (43) 

77/ 



^yit-' (44) 



leading to 



This is a factor of order 7r//3~^ times the simple estimate Eq. (30). Taking for definiteness 
1 < 7 < 3 and 1 < /? < 3 leads to 



^^a = (9to 800)(^^^^) A-^ (46) 



1.18 



In other words, taking all of the uncertainties into account Davis' proposal increases the 
simple estimate by a factor of order 10 to 1000. 

The proposal of Harari and Sikivie |3^, 35] is that the emission spectrum is flat, 

VemUt,k)=f{t) {^<^<f,) (47) 

a a 



This gives 



r.{t,k) = V/aV^lnf— ^) (48) 



n = ..yit-^r 1„(^)^ (49) 

= h^t-' (50) 

This is only a factor or order 7 times the simple estimate. Taking 1 < 7 < 3 it corresponds 
to 

/ f \^-^^ 
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In the two right hand columns of Figure 1 are shown the bands which correspond to these 
estimates. The fuh hnes correspond to T^eh = 10^^ GeV, the biggest possible value. The 
dotted lines (and their straight continuation) correspond to the absence of supersymmetry, 
or equivalently to Treh < 10^ GeV. The allowed window for ma is shown in the two left 
hand columns of Figure 2 as a function of Treh- 



Critical assessment of the proposals 

According to both estimates, most of the axions have angular frequency not far above ^^min 
(Eqs. ( [4^ and (^)). In other words the condition that their wavelength is much less 
than the string spacing, necessary for the axion concept to make sense, is not terribly well 
satisfied. Indeed, if 7 = 1 and (in the Davis case) /? = 1, it is at best marginally satisfied; 
most axions then have wavelength of order the string spacing. In that case, one might ask 



why the rough estimate Eq. (30) is not reproduced, which after all started with precisely 



the assumption that the wavelength was of order the string spacing. A clue to the answer 



appears when one calculates the mean square axion field [26|, 



= V.m{^) f (53) 

"-mm ^ ' 

(which follows from the fact that a relativistic plane wave with amplitude 9 and frequency 
uj has energy density f^O^u)"^). The corresponding mean square amplitude of the oscillation 
in 9 is 

02 = ^ ~ (.40)2 (Harari-Sikivie) (55) 
27r 

If we set 02 = 02 ^ pg[ ^ then Eqs. ( |45| ) and ( |5TD are roughly reproduced for the 

the case /3 = 7 = 1 (to be precise the result is a factor S/Svr smaller). In other words, the 
essential reason why the Harari-Sikivie proposal reproduces the simple estimate is that it 
leads to an amplitude ~ 1, and the essential reason why the Davis proposal is bigger by 
a factor of order r] is that it leads to an amplitude ~ 77^/2. 

This last fact is rather worrying. If i] had been a few orders of magnitude bigger than 70 
the amplitude would have been much bigger than 2tt. Such an oscillation amplitude would 
not in itself be unphysical, but the problem would be that it cannot be generated by the 
smooth, only marginally relativistic string oscillations envisaged in Davis' proposal; rather. 



such oscillations generate a typical amplitude ~ 1 |34|. Thus Davis' proposal would by 
inconsistent with the scaling assumption if 77 were much bigger than its actual value. One 
wonders what mechanism resolves the conflict when it arises, and whether this mechanism 
is already beginning to operate in the marginal case that we are dealing with.0 

^The above viewpoint is somewhat difFerent from the one taken by the present author in jiij. There, 
a less careful calculation gave a somewhat larger amplitude. More importantly, the fact that Davis' string 
oscillations cannot generate an amplitude 6 ^ 2n was not appreciated, and an estimate was made of the 
axion density resulting from such an amplitude. While this estimate would indeed be valid if such an 
amplitude were somehow generated, it does not apply to the the case at hand. I am indebted to Rick Davis 
for a clarifying discussion on this question. 
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Nevertheless, the statement that the oscillating strings emit radiation with a wavelength 
of order their spacing seems rather natural and to this extent the proposal of Davis is perhaps 
a reasonable working hypothesis. By contrast the proposal of Harari and Sikivie is extremely 
radical, in that it postulates the existence of a mechanism whereby the strings emit radiation 
impartially with all wavelengths, right down to the string thickness fa~^ < 10^"'^'^ GeV^^, 
no matter how late the epoch. As the string thickness is only 10"'^^ of the string spacing by 
the time radiation finishes, one would like to know something about the mechanism before 
taking the proposal too seriously 

As only strings as opposed to domain walls are involved it would seem feasible to perform 
numerical simulations. They would hopefully confirm the scaling behaviour, and determine 
the string spacing 7 defined by Eq. (31) and emission spectrum, thus allowing one to 
calculate on the assumption that axion number is conserved after wall formation. We 
end this section by briefiy asking about the validity of that assumption. 



After domain wall formation 

The string-wall network must annihilate well before the present to avoid cosmological disas- 
ter, and in order to make this possible we have assumed that the vacua on the two sides of 
a wall are identical (Eq. (p2[)). How does annihilation proceed? Based on |3^], the standard 
assumption in the literature seems to be that gravitational radiation is the only significant 
process. The underlying idea, inspired by what happens for a loop of gauge string, is that 
the network has significant structure only on macroscopic scales, of order the Hubble dis- 
tance and bigger. However, the walls are made out of the axion field and have a thickness 
of order l/m-a, and it therefore seems clear that collisions between pieces of wall will result 
in the emission of marginally relativistic axions at some level. Since gravity is so weak, 
this process is presumably the dominant one. As discussed in |26] wall annihilation gives a 
contribution 

""-(¥) (io4v) (^«> 

where tann is the wall annihilation time. Thus, axions emitted by the walls are probably 
significant if the Harari-Sikivie hypothesis is correct, but may not be if the Davis hypothesis 
is correct. 

It seems difficult to go beyond these rather crude considerations except through numer- 
ical simulations, which after wall formation would become extremely difficult. 



4.3 The no-string scenario 

So far we have assumed that axionic strings are generated in the early universe, either 
after infiation by the Kibble mechanism or during inflation by the quantum fluctuation. 
The criterion for the former possibility is T^eh ^ fa, and from Figure 3 one sees that this 
condition is rather similar to the condition Treh ^ ?decoup which corresponds to the presence 
of thermal saxinos in the universe. The latter condition can hold without the former only if 
Tj-eh is rather low, which means that the Kibble mechanism will certainly generate strings 
if saxinos generate a significant amount of entropy. 

The criterion for the generation of strings by the quantum fluctuation is less certain. 
Considering only the usual non-supersymmetric Peccei-Quinn held and ignoring any inter- 
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action with the inflaton field or the spacetime curvature, the criterion was estimated in 
to be Hi > 1 X /a, where Hi is the Hubble parameter when the observable universe leaves 
the horizon. This estimate should be fairly robust as a rough criterion, and may perhaps 
be thought of as a quasi-thermal effect, due to the Hawking temperature Hi/2tt. It implies 
that string production does not occur if Hi < 10^'' GeV (because /a > lO^'^GeV). While 
such a low value of Hi is possible [p^, 16|, most models of inflation require Hi ~ 10^3 to 



10^^ GeV, in which case the no-string scenario requires fa > 10^'^ GeV (Figure 3). 

A rather complete treatment of the no-string scenario has been given in and we 
recall briefly the main results, again ignoring any coupling of the Peccei-Quinn field to other 
fields. The axion field is homogeneous before the epoch T, except for an inhomogeneity 
which originates as a vacuum fluctuation during inflation. Writing 9 = 6 + 69, where 9 is 
the average of 9 in the observable universe, Eq. (^) for l^a becomes 

/ f \ 1-18 



n^ = 0.9xl0^-'^-^^j ^-\9' + <yt) (57) 
where cig, the average of (69)^, is given by 

''-^71 '''' 

The inhomogeneity of the axion field causes a primeval isocurvature density perturbation, 
with a flat spectrum given by [p^] 

^ 16 (^2 + ^2)2 (59) 



(It might cause bags of domain wall as well [23|, but this possibility has yet to be fully 
explored.) If present, a primeval isocurvature density perturbation contributes to the large 
scale microwave background anisotropy. Its effect on the rms anisotropy is the same as that 
of an adiabatic density perturbation whose spectrum 6jj at horizon entry is given by 

5h = 6x (2/15)^^/' (60) 
The COBE measurement implies |16| that 6h = (1-7 =b 0.3) x 10"^ which corresponds tc0 



Vli^ < 1.7 X 10-5 (61) 

At fixed 9, Eqs. (|^), (|^) and ( |6l| ) forbid a substantial region of the rrisrHi plane p8| . 
The forbidden region is minimised for ^ = 0, but there is no reason to expect 9 to be very 
small, unless a very small value is demanded by the requirement f^a < 1- Suppose that one 
sets 9 = min(.l7r, 9^) where 6q, is the value of 9 which makes i^a = 1, and fixes Hi. Then an 
interval of ma is forbidden, which covers the entire regime Hi < f^, < nipi'ii Hi > 10^^ GeV 
but which shrinks to nothing for Hi < 10^'^ GeV. 

The right hand column of Figure 2 illustrates the case Hi < 10^^ GeV; to the left of the 
diagonal line there is no constraint on rria except the astrophysical one. As Hi is increased 
from lO^'^ GeV to 10^^ GeV a forbidden region rapidly develops and the situation becomes 
practically the same as for the string scenario. 

weaker constraint V^^^ < 10~* was used, corresponding to the use of pre-COBE quadrupole data, 
but the difference is not very significant for the present purpose. 
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5 Conclusions 



The calculation presented here is valid under the assumption that the correct extension of 
the Standard Model respects Peccei-Quinn symmetry and supersymmetry, the latter being 



implemented without introducing dramatic new features such as a very light gravitino |13|. 
Subject to this requirement, the calculation shows that saxinos of thermal origin can dilute 
estimates of cosmological mass densities by a factor of up to 10^, and non-thermal saxinos 
might cause even more dilution. Let us close by briefly considering the implication of a 
dilution factor for some commonly discussed cases [j^], leaving aside the axion which has 
already been treated. 

The baryon density is usually considered these days to be generated at the electroweak 
transition, and at least the more conservative proposals do not generate significantly more 
than the observed density (eg. |3^). Accordingly these models could not tolerate a large 
dilution factor, but as the theoretical situation is still very fluid it is probably too early to 
say anything very definite. 

Apart from the axion a favourite dark matter candidate is a massive neutrino species. 
Species with mass mu < 10 MeV decouple while still relativistic, and without any dilution 
factor 17 < 1 requires nii, < 100 eV. Decoupling occurs just before nucleosynthesis, so this 
bound is not affected by entropy production (which must finish by that epoch). Species 
with mass > 10 MeV on the other hand decouple while non-relativistic, and without 
dilution this opens up an additional allowed interval rriu > 1 GeV. In that case decoupling 
occurs at 70 MeV (mj^ /I GeV) [|l| so dilution can occur and it it does the allowed interval 
extends to lower masses. From Fig. 5.2 of [|l| one learns that with a dilution factor A 10^ 
the interval extends right down to rrii, ~ 10 MeV. 

Finally, consider a particle species with much weaker interactions than the neutrino, so 
that decoupling occurs at T ^ 1 GeV when g^, ~ 100. Without entropy dilution the mass 
needed to make = 1 is m ~ 1 keV, but with dilution the mass is increased by a factor A 
and so could be of order 10 MeV. Such a particle species would constitute cold dark matter, 
instead of 'warm' dark matter as in the usual case. 

From these examples one sees that a significant entropy dilution factor could have pro- 
found implications. Whether or not such a factor is generated by thermal saxinos depends 
mainly on whether the reheat temperature is high or low, and we are reminded that an 
understanding of this quantity is one of the outstanding problems of theoretical cosmology. 
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Figure Captions 

FIG. 1. — Plots against axion mass. All scales are logarithims to base 10, and in each 
plot the horizontal axis is ma/leV. The rows correspond to different saxino masses as 
indicated. The left hand column shows various temperatures in units of 1 MeV, the full 
line indicating the temperature T after which axion number is conserved, and the dotted 
lines indicating the three temperatures Tdecay < r= < Tsaxeq defined in the text (entropy 
generation occurs only in the narrow interval Tjecay < T < T=). The second column gives 
the entropy production factor A, given by Eq. ( |T7|) or (dotted line) Eq. (|l6|). The two right 
hand columns give the axion density f^a, for T^eh = 10^^ GeV (full lines) and Tj-eh < 10^ GeV 
(dotted lines and their continuation). In each case a band of values is shown to take into 
account the uncertainties discussed in the text, and the second case is identical with the 
no-supersymmetry result. The hatched regions correspond to the nucleosynthesis constraint 
and (for the two right hand columns) the constraint Qg, < 1. 
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FIG. 2. — The axion window. In each plot the vertical axis is log]^o("^a/l eV) and the 
horizontal axis is log2o(^reh/l GeV). The first row corresponds to the non-existence of 
the saxino and the others to different saxino masses. The three columns correspond to 
different possibilities described in the text. The upper cross hatched region is forbidden by 
accelerator physics and astrophysics. The lower cross hatched region is forbidden either by 
r^a < 1 or by the requirement that saxino decay must not interfere with nucleosynthesis. 
The latter requirement is relevant only for the second and third rows, giving the right hand 
horizontal part of the boundary in each plot. In each plot (except on the first row) the 
diagonal full line and its straight continuation is the line below which thermal saxinos are 
absent, so that there is no dilution factor (non-thermal saxinos are not considered in the 
present paper). The diagonal dotted line in the right hand column is the line below which 
the Kibble mechanism does not produce axionic strings. 

FIG. 3. — Various regimes of the ma-Treh plane, as discussed in the text. The horizontal 
line marks the regime Hi > /a for the case Hi = 10^^ GeV, and it moves up one decade for 
every decade that Hi moves down. 
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